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1 Introduction 

There has been recently an increasing interest in infrared modifications of gravity which 
could, in principle, explain the accelerated expansion of the Universe without introducing 
a "dark energy" component. However, before these models are taken seriously, their con- 
sistency has to be checked both from theoretical point of view and by comparison with 
the experimental data. One class of such models involves spontaneous breaking of Lorentz 
symmetry by the space-time dependent condensates of scalar fields coupled to gravity 
through derivative couplings [1, 2, 3]. It has been shown that such models may have a non- 
pathological behavior in the perturbative regime and may exhibit infrared modifications 
of the gravitational interactions. The models of this class are naturally called "massive 
gravity" since the graviton generically acquires a non-zero mass due to interactions with 
the scalar fields [3, 4]. An important feature of these models is that they are formulated 
in a non-perturbative way which makes it possible to study non-linear solutions such as 
those describing modified cosmological evolution [5] or black holes [6, 7] (for a review and 
further references see, e.g., Ref. [8]). 

In the conventional General Relativity (GR), the existence of solutions with horizons 
(e.g., black holes) raises questions of consistency of a more general kind, namely, the con- 
sistency with the general laws of thermodynamics. It has been argued [9, 10, 11] that in 
GR such consistency can be obtained by assigning black holes a certain temperature and 
entropy. The temperature of a black hole characterizes its Hawking radiation, so that the 
black hole would be in a thermal equilibrium with the thermal bath at this temperature. 
The entropy of a black hole was argued to be proportional to its horizon area. In this way, 
the net entropy of the black hole and the outer region never decreases, generalizing the 
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second law of thermodynamics to processes that include black holes. These thermodynam- 
ical properties of black holes in GR are believed to be connected to fundamental principles 
of quantum physics, such as unitarity. The validity of the thermodynamical description of 
black holes should then give us some insight into quantum aspects of gravity. For example, 
in string theory for a certain class of extremal black holes the Bekenstein entropy has been 
reproduced by counting the microscopic states of the compact objects [12, 13]. 

A question that arises naturally is how the thermodynamic properties of the black holes 
are changed in the modified gravity models, in particular, in massive gravity, and whether 
these changes preserve consistency with the thermodynamic laws. In the context of the 
ghost condensate models [2] this question has been raised in Refs. [14, 15]. In Ref. [14] a 
gedanken experiment involving a black hole has been proposed which allows the transfer of 
heat from a cold body to a hot one, thus violating the second law of thermodynamics (for 
the discussion of the subtleties of the arguments see refs. [16, 17]). In Ref. [15] the violation 
of the second law of thermodynamics was related to the presence of negative energy states. 
On the contrary, in the context of the TeVeS models [18] it was conjectured [19] that the 
second law of thermodynamics holds if the effective graviton radiation temperature and 
the Hawking radiation temperature are equal. 

The aim of this paper is to check the validity of the laws of thermodynamics in massive 
gravity by making use of the exact black hole solutions of Ref. [7] (solutions of this type 
were first found in bi- metric models in Ref. [20]). Unlike conventional black holes, these 
solutions depend on two parameters: the mass M and the "scalar charge" Q. We will 
concentrate on solutions which can be interpreted as (modified) black holes, i.e., have a 
finite ADM mass (and, therefore, produce the Newtonian 1/r potential at large distances) 
and have an event horizon. We will first find the expressions for the temperature and 
entropy of the modified black holes. This will be done by calculating the surface gravity 
and by making use of the Wald's representation of entropy as a Noether charge. Having 
obtained the explicit expressions for the entropy and temperature in terms of the black hole 
mass and scalar charge, we then consider the process of merging of two black holes. As we 
will argue, in the limit when one of the black holes is much smaller than the other, the mass 
and the scalar charge of the resulting black hole is the sum of masses and charges of the 
constituents. Thus, we will be able to compare the entropies of the initial and final states 
and check the second law of thermodynamics. We will see that when the negative mass 
states are present, one can decrease the total entropy during the merger. Moreover, one 
can create a zero-temperature black hole — an extremal state which, by an infinitesimal 
change of initial parameters, can be converted into a state with naked singularity. Thus, 
the second and third laws of thermodynamics can be violated. 

The paper is organized as follows. In Sect. 2 we review the static spherically symmetric 
solutions of massive gravity. In Sect. 3 we first compute the temperature and entropy of 
modified black holes in massive gravity. We then verify the second and the third laws of 
thermodynamics and show that both can be violated. Sect. 4 contains the summary of the 
results and their discussion. 
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2 Static Spherically Symmetric Solutions in Massive Gravity 



The massive gravity model which is used in this paper is described by the following action 
[3]: 

S = j dx^./^ [-Ml{R + A^J^] (2.1) 

where the first term is the usual Einstein-Hilbert action and the second one is a certain 
function (to be specified below) of the space-time derivatives of the four scalar fields 

minimally coupled to gravity. This model should be viewed as a low-energy effective 
theory with the cutoff scale A. 

The model (2.1) generically admits a flat-space "vacuum" solution 

The vacuum possesses rotational symmetry provided the function T is invariant under the 
rotations of the fields in the internal space. The Lorentz symmetry is, in general, broken. 
Requiring that the action (2.1) is invariant under the following symmetry, 

</.^^</.^ + H^(/), (2.2) 

where H* are arbitrary functions of (fP , ensures that perturbations about this solution 
contain only two propagating degrees of freedom [3] which are two polarizations of a massive 
graviton with the mass of order m = A^/Mpj, where Mpi is the Planck mass. The model 
does not contradict the most obvious experimental constraints [5, 21, 22] for graviton 
masses as large as 10"^'' eV. 

The ansatz for the static spherically symmetric solution can be written in the following 
form [7]: 

ds"^ = a{r)dt^ - /3(r)dr2 - (d^^ ^ gij^2 ^^^2^ ^ 
/ = A2 [t + h{r)], 

ct>' = ct>{r) , (2.3) 

r 

where the coordinate transformations r ^ r' = r'{r) and t — )■ t' = t + T(r) have been used 
to eliminate some of the fields. 

The analytical black hole solutions can be obtained if the function T has the form 

12 /2 \ o 
J^=^ + 6(- + lju;i-t/;? + ?,wiW2 - 2w^ + 12, 

where A is a positive constant and 

Note that the dependence of J- on the scalar fields through two combinations X and W^^ 
ensures the symmetry (2.2). 
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Figure 1. Newtonian potential 2$ = goo — 1 in the three different cases: M > and Q > 0, M > 
and Q < 0, M < and Q > 0. For eomparison, the potential of the Schwarzschild solution of the 
mass \M\ is shown by the dashed line. 



With this function J^, the black hole solution reads 

2MGn Q 



a{r) = 1 
I3{r) 



1 



h{r) 



a{r) 

± I ^ 

a 



^ ^, I Q A(A - 1) ^ ^ 



1/2 



r = r. 



(2.4) 



Here M and Q are two arbitrary integration constants. We will assume in what follows 
that A > 1. In this case the asymptotics of the gravitational potential is Newtonian and is 
given by the parameter M which determines the ADM mass of the solution. 

The solution may possess an horizon which, in this case, is given by the largest root 
of the equation 

2MGn , Q 



(2.5) 



Depending on the signs and relative values of the parameters M and Q, there are three 
different cases when the horizon exists: (M > 0, Q > 0), (Af > 0, Q < 0) and (M < 0, Q > 
0).The three corresponding solutions are shown in Fig. 1 

In the case M > and Q > 0, the black hole has an attractive gravitational potential 
at all distances. The attraction is stronger than that of the usual Schwarzschild black hole 
of mass M; The horizon size of the modified black hole is also larger (left panel of Fig. 1). 

When M > and Q < 0, the horizon only exists when the condition 



2MGn > \ \Q\ 



l/A 



1 



A-1 



A-l 
A 



(2.6) 



is fulfilled. In this case, the Newton's potential is also always attractive. However, the 
attraction is weaker than in the Schwarzschild case, and the horizon size is smaller. 

Finally, in the case M < and Q > the Newton's potential is repulsive at large 
distances and attractive near the horizon. This case is interesting since it doesn't exist in 
GR. Although formally the black hole mass in GR can also be taken negative, this does 
not correspond to a physical solution because of the naked singularity. Another reason 
to disregard such solutions is the null energy condition which holds for the matter stress 
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tensor [23, 24, 25, 26]. Neither of the arguments exist in the modified gravity case. At 
small distances the repulsion changes to the attraction, which creates the event horizon 
hiding the singularity (right panel of Fig. 1). Also, massive gravity doesn't satisfy the null 
energy condition, allowing for negative mass states to be constructed, e.g., as in the ghost 
condensate model [27]. 

3 The Thermodynamical Properties of the Black Hole Solutions 

In order to analyse the validity of the second law of black hole thermodynamics we first 
compute the temperature and entropy of the modified black holes. We use the approach 
by Wald [28, 29, 30, 31, 32] based on the Noether's charge. As we will show, the entropy 
is given by the Bekenstein-Hawking formula, that is equal, in the appropriate units, to the 
one quarter of the horizon area. 

3.1 The Black Hole Temperature 

Due to quantum particle creation, black holes emit thermal radiation [11]. The temperature 
of this radiation is determined by the surface gravity k of the black hole, that is, the 
acceleration experienced by a test body at the black hole horizon. The event horizon of 
a static spherically symmetric black hole corresponds to a killing horizon, i.e, a surface to 
which a killing vector field is normal (bifurcation surface). The surface gravity k, at any 
point of a killing horizon % is defined by 

e^ae = <\ (3.1) 

where is the killing field normal to %. For a static spherically symmetric metric of the 
form 

ds^ = a{r)dt^ ^dr^ - (d^^ + sin^ Odcf) , (3.2) 

a{r) 



the surface gravity equals 

"'^^^ (3.3) 



r=rH 

where rn is the horizon radius. The zeroth law of black hole mechanics asserts that k 
is constant over the event horizon of a stationary black hole. Although in GR one needs 
to use the Einstein's field equations to prove this statement, and thus its generalizations 
to other theories of gravity is questionable, the zeroth law trivially holds for spherically 
symmetric black holes. 

Making use of the explicit solution (2.4) one finds 

K 1 1 X Q 



2tt Airrn 47r r^+^ 

1 / 1 + (A - l)Qr^n ^ 
4vr \ rn ' 



(3.4) 



One recovers the temperature of the Schwarzschild black hole in the limit of zero scalar 
charge, in which case Th = l/AiirH as expected. Moreover, since the temperature of the 
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f max 



Figure 2. The red line corresponds to the temperature of a black hole with negative scalar charge. 
It reaches its maximum at rmax = [(-'^^ ~ 1)IQI] The blue line represents the temperature 
as a function of the horizon radius for positive scalar charge solutions. The dashed line is the 
Schwarzschild case. 

black hole corresponds, from a mathematical point of view, to the tangent of the Newton's 
potential at the event horizon, it is easy to conclude from Fig.l that the existence of an 
event horizon implies Th > 0. 

Interestingly, the Hawking temperature behaves differently for positive and negative 
scalar charges. For positive Q, the temperature is larger than at Q = and decreases with 
rn as in the case of a conventional Schwarzschild black hole. On the contrary, at Q < 
the temperature is smaller than in the Schwarzschild case. Moreover, its dependence on 
is not monotonic, with a maximum reached for rn = ?'max = [(^^ — ^^'^ (see Fig. 2). 

When rn > r^a.x, the specific heat is negative, as in the Schwarzschild case, corresponding 
to the fact that the black hole heats up as it radiates energy, while when rn < ?^max) 
the specific heat is positive, which means that the black hole cools down as it radiates, 
eventually reaching an equilibrium. 

Another interesting observation is that the temperature reaches zero at a finite event 
horizon radius for Q < 0. This is closely related to the third law of black hole thermody- 
namics and will be discussed in the next section. 

3.2 The Noether's Charge as the Entropy 

As has been proven in Refs. [28, 29, 30, 31, 32] by assuming a theory which admits black 
holes with a bifurcate horizon and a well-defined mass at infinity, the entropy of the black 
hole is given by 

Sh = 27r [ Q[^], (3.5) 

where Q[i^] is the Noether charge related to the Killing field normal to the horizon T-L. 
For a general theory with matter fields "0 and the Lagrangian 

C= C {gah, TZabcd, -0, V', ^ {ai,...^ a,)i^^ , (3-6) 
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this quantity can be recast in the fohowing form [29]: 



Sh = -27T [ 



d 



,D-2 




5C 



(3.7) 



JH 



where a is the determinant of the metric on Ti and iab is the binormal vector to the 
bifurcation surface. The black hole entropy is then related to the Noether's charge of 
diffeomorphism under the Killing vector field which produces the event horizon. 

Eq. (3.7) can be directly applied to the case of massive gravity with the action (2.1). 
For metric tensors of the form (3.2), the relevant Killing vector is dt, while the binormal 
vector iab has the following components: etr = —£rt = 1 and the other components vanish. 
The explicit expression for the entropy is then 



where the factor 4 is a consequence of the antisymmetry property of the Riemann tensor 
and the binormal vectors. The index (s) is to emphasize that the functional is evaluated on 
the solution (2.4). Substituting the Lagrangian (2.1) and making use of the fact that the 
function does not dependent on the Riemann tensor (the metric and the Riemann tensor 
are treated as independent variables in eq (3.8)), we arrive at the Hawking-Bekenstein 
formula : 



where Ah is the area of the black hole horizon in the Planck units = 1. 

As it has been shown in Ref.[28], this expression for the entropy is automatically 
consistent with the first law of thermodynamics. It can also be checked directly by making 
use of the expressions (3.9) and (2.5) at dQ = 0. 

3.3 The Second and Third Laws of Thermodynamics 

Let's first examine the generalized second law of thermodynamics. To this end consider 
the change of entropy in the process of coalescence of two black holes characterized by the 
masses Mi and m2 and scalar charges Qi and q2, respectively. Let both black holes be 
large enough so that the Hawking radiation can be neglected and their horizon radii are 
larger than the inverse cutoff scale, in which case the solutions are within the region of 
validity of the effective theory. 

Before the coalescence, when the interaction of the two black holes can be neglected, 
their entropy is simply the sum of the entropies of the two isolated black holes as given 
by eqs. (3.9) and (2.5). In order to check the second law of thermodynamics we need to 
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determine the entropy of the final state. Let us argue that in the limit when one of the 
black holes is much larger than the other, and the scalar charges of the black holes are not 
parametrically larger than their masses, 

|Mi|»|m2|, |Qi|'/^< |Mi|, k2|'/^<|m2|, (3.10) 

the result of the coalescence is a black hole with the mass Mi + m2 and the scalar charge 
Qi + q2- Then the entropy of the final state is given by eqs. (3.9) and (2.5), and the net 
change of entropy can be easily calculated. 

The argument consists of two parts. First, let us show that in the above limit the 
total energy of gravity waves emitted in the process of coalescence is parametrically small 
as compared to the mass of the smallest black hole. Assume that the characteristic time 
scale of the coalescence is much smaller than the inverse graviton mass (this can always 
be arranged in view of the tiny value of the latter). Then the gravitational waves can be 
considered massless. The metric perturbation hij in the case of the quadrupole radiation 
is estimated as follows: 



where Qij is the second time derivative of the quadrupole moment and r is the characteristic 
size of the system. The energy density in the gravitational waves is, therefore, of order 



-XX- (3-12) 
pi 



where u is the characteristic frequency of the emitted waves as seen by an asymptotic 
observer. Setting the size of the system to r ~ Rh and the frequency to a; ~ 1/Rh, Rh 
being the horizon size of the large black hole, we have 

Q ~ ijJ^Q ~ uP'm2R?H 

and, therefore, the total energy emitted over the time period ~ Rh is 

^.a. - -2 (^) ~ -2 (^) , (3.13) 

where rn is the horizon size of the small black hole and we have assumed that the presence 
of scalar charges does not change completely the horizon sizes of the two original black 
holes. We see from eq. (3.13) that in the limit m2 ^ Mi the gravitational radiation during 
the coalescence of the two black holes can be neglected. 

We now turn to determining the mass and the scalar charge of the resulting black 
hole. This can be done by considering the asymptotic gravitational potential created by 
the coalescing black holes and extracting the coefficients in front of the 1/r and 1/r^ 
terms. These coefficients determine the mass and the scalar charge of the resulting black 
hole, respectively. 

To clarify the logic, consider first the coalescence of the two black holes of masses Mi 
and 777-2 in the conventional GR. In this case the asymptotic gravitational potential of each 
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black hole is determined by its mass and satisfies the linear equation. Thus, before the 
coalescence and to the leading order in 1/r, the gravitational potential of the two black 
holes is given by the sum of the potentials of the individual black holes. The coefficient 
of 1 jr is independent of the distance between the black holes (as long as this distance is 
much smaller than r) and is given by the sum of the black hole masses. Since after the 
coalescence this coefficient is determined by the mass M of the resulting black hole and 
there is no substantial emission of the gravitational waves, we conclude that M = Mi + m2- 
In the case of the modified black holes the situation is more complicated because the 
corresponding solution is non-linear at all distances. This introduces a correction to the 
simple addition of masses and charges, 

M = Mi+m2+5m- Q = Qi+q2+ Sq. 

However, one can choose the parameters Mi, m2, Qi and q2 in such a way that these 
corrections are negligible. Making use of an exactly treatable spherically symmetric case 
to estimate 6q and 5m (see Appendix A), one may argue that this can be achieved by 
requiring that 

• (3-14) 



The last condition can always be satisfied for a sufficiently small value of the graviton mass 
H without leaving the region of validity of the effective theory. In this case, the effects of 
non-linearity on the asymptotic potential can be neglected and the arguments presented 
above imply that the resulting black hole has, to a good accuracy, the mass M = Mi + m2 
and the scalar charge Q = Qi + q2- 

Having determined the mass and the scalar charge of the final black hole, we can now 
calculate the change of the entropy in the process of coalescence, 

AS'feh = S'final — S'initial = ['IZjj — Rjj - rjj) 

^7T{nl-Rl), (3.15) 

where TZh is the radius of the event horizon of the final black hole. 

In the limit (3.10) the entropy change can be calculated expanding in powers of m2/M\. 
To the leading order, the answer is linear in the mass of the small black hole, 

^ '""^^"^^ , = ^, (3.16) 

1 + (A - l)Qii?^^ Th' ^ ^ 

where Th is the temperature of the large black hole. Thus, the entropy may increase or 
decrease, depending on the sign of the mass of the small black hole. For negative masses, 
the generalized second law of thermodynamics is violated. 

The decrease of the entropy when negative mass black holes are involved can be checked 
explicitely in the particular case of A = 2. In this case eq. (2.5) for the horizon in terms of 
the mass and the scalar charge can be solved exactly. 



Rh = M + ^/aP + Q. (3.17) 
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The existence of the horizon requires + Q > 0, which is always the case for a positive 
scalar charge. For negative Q, the condition (2.6) has to be imposed. 

The change of the black hole entropy can be computed explicitly giving 



ASbh =2tt 



(Ml + m2)\/(Mi +m2)2 +Qi + 



92 



+ 2Mim2 - MiJMf + Qi- miJm^ + q2 



(3.18) 



One can see from this expression that the entropy always increases when both black holes 
have positive masses. However, when at least one of the masses is negative and Mi ^ 
\m2\, one arrives at the same conclusion as before: there is a decrease of entropy and the 
generalized second law of thermodynamics is violated. 

Consider now the third law of thermodynamics. In application to black holes, it states 
that a black hole with a non-vanishing temperature cannot reach zero temperature in a 
finite sequence of operations [33]. In our case, the zero temperature corresponds to an 
"extremal" situation when the inequality in eq. (2.6) becomes an equality. Indeed, the 
black hole temperature can be expressed in terms of the surface gravity which, according 
to eq. (3.3), is proportional to the slope of a{r) at the horizon. It is then clear from Fig. 1 
that the zero temperature can only be reached in the case M > 0, Q < (represented on 
the middle panel) when the two roots of eq. (2.5) coincide. Thus, an extremal black hole 
satisfies the condition 

2M = -f\Qf (3.19) 



with j3 = 1/X and 



A-l 
A 



Note that both 7 and /3 are positive real numbers. For simplicity, we will concentrate on 
the case A = 2 from now on. In this case, (3 = 1/2 and 7 = 2. The generalization to other 
cases is straightforward. 

Consider again the coalescence of two black holes. Let the large black hole be nearly 
extremal, so that Qi < and 

Ml = {l + e)\Qi\^l\ 

where e is a small positive number (e = corresponds to the extremal black hole). If the 
small black hole has a positive mass and a scalar charge satisfying the condition (2.6), it 
is easy to see that the black hole resulting from coalescence will also satisfy the condition 
(2.6) and thus will have a non-zero temperature. Indeed, the final state is characterized by 
the inequality 

(Ml + + Qi+q2 = {Ml + Qi) + {ml + ^2) + 2m2Mi > 0. 

On the contrary, if the small black hole has a negative mass, one can always choose its 
charge such that it converts the large black hole into an extremal one. For example, the 
mass of the small black hole can be taken m2 = — 2e|Qi|^/^. The condition m2 ^ Mi is 
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still fulfilled for small enough values of e. The scalar charge of the second black hole can 
be chosen as q2 = 5\Qi\ with 5 a small positive number. Then, the final state is a black 
hole with mass M. = (1 — e)|Qi|"^/^ and a scalar charge \Q\ = (1 — (5)|(5i|, which means that 
= |Q|^/^ if we choose (1 — e) = (1 — dY^"^ . A non-extremal black hole has turned into 
an extremal one in a one-step process. Thus, when negative mass black holes are present, 
the third law of thermodynamics is not fulfilled. 

The third law of black hole thermodynamics is closely related to the cosmic censorship 
conjecture [34], which states that all singularities resulting from a gravitational collapse 
are always hidden by the horizon. With an obvious modification, the above process that 
leads to an extreme black hole can be used to create a naked singularity. This possibility 
is again related to the existence of the negative-mass black holes and, ultimately, to the 
violation of the null energy condition in massive gravity. 

4 Discussion 

In this paper we have addressed the validity of the generalized laws of thermodynamics 
in massive gravity models. In these models the black hole solutions are modified by the 
presence of the scalar "hair" . The analog of the Schwarzschild black hole - the spherically 
symmetric solution - depends on two parameters, the mass and the "scalar charge" char- 
acterizing the hair strength. The presence of two free parameters makes the asymptotics 
of the gravitational field of the black hole essentially independent of its behavior near the 
horizon, allowing for negative-mass solutions without naked singularity. 

Making use of the exact black hole solutions, we have constructed explicit examples 
of the black hole mergers which violate the generalized second and third laws of thermo- 
dynamics. The existence of such processes is in accord with the general theorems of the 
black hole thermodynamics which require the weak energy condition to hold. The latter 
condition, together with the null energy condition, is violated in massive gravity models. 
Indeed, the examples we have constructed involve the negative-mass black holes. Note that 
the violation of the null energy condition in massive gravity is related to the presence of 
the superluminal modes [35] which are likely to be responsible for the existence of the black 
hole hair [6]. 

Although the negative mass black holes are classical solutions of massive gravity which 
exist on the same footing as the conventional positive-mass black holes, this does not 
guarantee that they can be created from positive-mass states in the course of the evolution. 
If this were the case, the violation of the second law would imply that one may devise a 
process involving black holes which would convert heat into work. Moreover, the violation 
of the third law would mean that the cosmic censorship conjecture is not true, so that one 
can create naked singularities starting from singularity- free states. 

It is not inconceivable that the creation of the negative-mass states can be forbidden 
dynamically. In fact, it has been conjectured in the context of the ghost condensate model 
that an average null energy condition (ANEC) may prevent the entropy of the black hole 
from decreasing, in a coarse-grained sense [17]. This means that the second law of ther- 
modynamics may be violated locally, while in average (during long time periods) it holds. 
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The ANEC can be traced back to the presence of the conserved Noether's charge related 
to the shift symmetry of the ghost condensate field. In our case, a similar internal shift 
symmetry is also present, so one may wonder whether this symmetry implies ANEC which 
may prevent the breakdown of thermodynamical laws. 

The situation, however, is not exactly the same in massive gravity models. First, 
massive gravity possesses four scalar fields (fP and 0*, the former being an analog of the 
ghost condensate field. Even though these fields decouple from the ordinary matter, they 
are coupled to each other, so the scalar sector is more complicated than in the ghost 
condensate model. For this reason we have not been able to directly generalize the argument 
based on ANEC, except for linear perturbations above the Minkowski background. This, 
however, is not sufficient to argue against the formation of the negative-mass black holes, 
so the question remains open. 

Moreover, the original argument [17] is partially based on the fact that a negative 
Noether's charge is strongly disfavored since it is leading to UV instabilities. However, 
the massive gravity is free from UV instabilities, at least for backgrounds which can be 
approximated as flat in the UV limit (see [5] for details). This implies that part of the 
argument in [17] doesn't go through in our case. 

The final observation is that ANEC cannot protect the third law of thermodynamics 
and the cosmic censorship conjecture. Indeed, unlike the entropy which may decrease lo- 
cally but still increase in average, the naked singularity, once created, cannot be undone 
no matter what happens in other places. The breakdown of the cosmic censorship con- 
jecture may therefore occur even if ANEC holds. The creation of naked singularities may 
signal the problem of the model, unless they are somehow regularized by an appropriate 
UV-completion. 
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A Estimate of charge and mass correction 

Since the equations of massive gravity do not linearize asymptotically because of the slow 
decay of the field cfP [7], the state of two close black holes with parameters Mi, Qi and m2, 
(72 has asymptotic mass and scalar charge (i.e., the coefficients in front of 1/r and 1/r^ in 
the Newtonian potential) equal to 

M = Mi+m2 + 6m, Q = Qi + q2 + Sq. 

Our goal here is to estimate the corrections 5m and 6q and find the parameters for which 
these corrections can be neglected. 
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For the estimate we make use of the fact that for a spherically symmetric matter 
distribution the equations can be solved exactly. We model the large black hole by a sphere 
of constant density of the total mass Mi and radius close to the horizon size R ~ 2GnMi, 
and a small black hole by a spherical shell of the same density and the mass m2, covering the 
sphere R. Even though the geometry is different, we expect that our model configuration 
reproduces correctly the parametrical dependence of 6m and 6q on R, Mi and m2. 

In a spherically symmetric case it is straightforward to obtain the scalar charge as of 
a sphere of a constant density as a function of the mass Mi and radius R [36]. Matching 
boundary conditions, we obtain the following relation for the scalar charge 

Qi = C{\)MiGnIi^R^^^. (A.l) 

where C{\) is a constant depending on A which is of order one for A varying between one 
and two, and ^ is the graviton mass. 

Making use of eq. (A.l), one may compute the change of the scalar charge when a thin 
shell around the constant density sphere is added. The result is 

This estimate can be rewritten in terms of the parameters of the two black holes as follows: 

Requiring 6q <^ q2 we arrive at the condition (3.14) for the worst case scenario, i.e. when 
m2 ~ A. If the condition (3.14) is satisfied, it automatically implies that 5m <C m2- 
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